Landscape Predictions from Cosmological Vacuum Selection 
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In BP models with hundreds of fluxes, we compute the effects of cosmological dynamics on the 
probability distribution of landscape vacua. Starting from generic initial conditions, we find that 
most fluxes are dynamically driven into a different and much narrower range of values than expected 
from landscape statistics alone. Hence, cosmological evolution will access only a tiny fraction of 
the vacua with small cosmological constant. This leads to a host of sharp predictions. Unlike 
other approaches to eternal inflation, the holographic measure employed here does not lead to 
"staggering" , an excessive spread of probabilities that would doom the string landscape as a solution 
to the cosmological constant problem. 



I. INTRODUCTION 

The task of making predictions in the landscape of 
string theory represents an enormous challenge. We must 
survey the metastable vacua in the theory and find the 
relative abundance of various low-energy properties. But 
vacua without observers will not be observed. The great 
variability of low energy physics in the string landscape 
demands novel methods for capturing such selection ef- 
fects. 

Interposed between these two tasks is a third problem, 
which will be the focus of this paper: vacuum selection 
by cosmological dynamics. It is not enough for a vacuum 
to exist in theory. It will be relevant only if it can be dy- 
namically realized from generic initial conditions, either 
as regions in spacetime or branches of a wavefunction. 

The former viewpoint insists on treating spacetime 
semiclassically on a global scale. This has no operational 
meaning, since no observer can see beyond his own hori- 
zon Moreover, it leads to problematic infinities and 
ambiguities that have plagued the definition of measures 
in eternal inflation; see, e.g., Refs. @, i, i, i, H, 0, [1 
for discussions and some recent proposals. Thus, the 
global viewpoint seems poised to suffer the same fate as 
the ether: a false convenience increasingly recognized as 
a burden. (A similar conclusion was reached earlier by 
studying quantum aspects of black holes [H, |T(| , and it 
seems natural indeed that it should extend to cosmology.) 

One of us recently proposed a "holographic" mea- 
sure Q that refers only to a sin gle causally connected 
region, or causal diamond. In Ref. [IlJ , it was shown that 
this measure not only overcomes problems that plagued 
anthropic predictions of the cosmological constant, but 
is able to maintain this success when specific anthropic 
conditions are traded for the much weaker assumption 
that observers require merely free energy. 
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Here we apply a different aspect of the proposal of 
Ref. 0]. We focus not on anthropic selection, but on 
cosmological selection. We shall find that dynamical ef- 
fects can significantly suppress or enhance the probabil- 
ity of observing a given vacuum, and they can interfere 
with anthropic selection in some models. However, unlike 
another recently proposed measure p|, the holographic 
measure does not lead to excessively uneven probabil- 
ity distributions that would render the string landscape 
ineffective at solving the cosmological constant problem. 

We will follow Ref. 12] in modelling the string land- 
scape as arising from a large number of possible four-form 
flux configurations. We mainly consider a model with 
250 fluxes, with fixed charges of order 1/30. This yields 
a large number of metastable vacua, of which 10 121 have 
vacuum energy comparable to the observed value. Of 
these, however, we find that only 10 80 are accessed by a 
typical worldline, starting from generic initial conditions. 

The selected ensemble is characterized by 250 probabil- 
ity distributions over the integers, one for each flux. The 
distributions differ — in many cases, drastically — from the 
distributions one would have obtained by simply restrict- 
ing the landscape to vacua with small cosmological con- 
stant. Thus, cosmological selection leads to thousands 
of distinct predictions. Many correspond to probabilities 
that are so close to or 1 that even a single conflicting 
observation would rule out the model. 

The specific predictions we obtain apply only to the toy 
model we study, but they do allow us to draw more gen- 
eral lessons. Most importantly, our results demonstrate 
that cosmological dynamics can be a powerful constraint. 
It reduces the effective size of the landscape drastically, 
leading to a large number of strong predictions. Another 
general lesson is that fast decays happen first. In natu- 
ral models, it takes hundreds or thousands of tunneling 
events to get from a generic initial vacuum to a vacuum 
with observers. Any metastable internal configuration 
that gives rise to a relatively fast decay channel but is 
unlikely to arise late in the decay chain will not be ob- 
served. 

In the real string landscape, the selection effects will 
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likely be as strong as in the BP models studied here, 
but they may be harder to characterize in terms of a 
set of nearly independent flux distributions. For exam- 
ple, charges will not be independent of fluxes, and their 
dynamical effects can become intertwined. 

The specific flux configuration associated with our own 
vacuum is unlikely to be observed directly in the near fu- 
ture. However, fluxes are related to coupling constants 
and other low-energy parameters. It will be important 
to understand such correlations in realistic and represen- 
tative sets of string compactifications. Then it will be 
possible to understand the impact of cosmological selec- 
tion effects on the distribution of low-energy parameters. 

a. Outline Sec. |TT] includes brief descriptions of the 
BP model, and of the holographic cosmological proba- 
bility measure we use to compute cosmological selection 
effects. We focus on the question of how this measure can 
be effectively applied to a model with ten to the hundreds 
of vacua. This would seem to require inverting a matrix 
of the corresponding rank. However, we develop statis- 
tical methods that are both adequate and very efficient. 
(Details of the techniques are given in Appendices 1X1 and 
IBT ) We also discuss our choice of initial conditions, and 
we show that our results will be the same for any initial 
probability distribution that depends only on A and does 
not strongly favor small or negative values of A. 

In Sec. IIII1 we present our results. For a representa- 
tive model, we contrast the distribution of fluxes, and 
the number of vacua with small cosmological constant, 
before and after cosmological selection effects are taken 
into account. We point out specific predictions that can 
be made with high confidence, focussing on features that 
would have been either absent or strikingly different if 
cosmological dynamics had been neglected. We consider 
other models to illustrate that the strength of cosmolog- 
ical selection effects varies and that, in special cases, a 
model that naively contains enough vacua may fail after 
cosmological selection is accounted for. 

In Sec. HVi we try to give a qualitative explanation of 
the features seen in the cosmological flux distribution. 
We identify them as imprints of the dynamics at differ- 
ent stages in the decay chain, when different types of 
transitions dominate. 

In Sec. El we compare our results to those obtained 
by Schwartz-Perlov and Vilenkin [l3j from a different 
proposal for the probability measure 5]. In toy-BP- 
models with fewer vacua, the alternative measure leads 
to a "staggered" probability distribution, which would 
leave the string landscape effectively underpopulated and 
unable to solve the cosmological constant problem. The 
measure of Ref. [f| is difficult to apply explicitly in models 
with a realistic number of vacua, and to extend beyond 
first order in perturbation theory. However, we show in 
Appendix [C] that it is equivalent to the holographic mea- 
sure with particular (and from our point of view, highly 
unnatural) initial conditions. With the help of this trick, 
we are able to argue that the staggering problem is likely 
to persist at higher orders when the measure of Ref. [5[ 



is applied to large BP models. 



II. MODEL AND MEASURE 



A. Probability measure, course graining, and 
Monte Carlo chains 



Consider a worldline starting in some metastable ini- 
tial vacuum. This vacuum will eventually decay, giving 
way to a new vacuum. The cascade continues until the 
worldline enters a terminal vacuum, which does not decay 
further. According to the proposal of Ref. Q , the prob- 
ability Pa assigned to vacuum A is the expected number 
of times the worldline will enter vacuum A on its way 
down to a terminal vacuum. 1 For a more general initial 
probability distribution P^' , one sums over the result 
obtained from each initial vacuum, weig hted by Pj 0) . 

The probabilities Pa form an A/"-dimensional vector, 
where Af is the number of vacua in the landscape. It 
satisfies the matrix equation Q 



(1 - T])P = r]P {0) 



(1) 



where t\ab is the relative decay rate, or branching ratio, 
from vacuum B into vacuum A. It can be obtained from 
the decay rate per unit time, kab, by normalizing each 
column of kab , 



Vab 



kab 



J2c K CB 



(2) 



except for columns corresponding to terminal vacua 
(which vanish both for n and for 7/). 

In a large landscape, one expects that both the ini- 
tial distribution, , and the distribution resulting from 
dynamical selection, Pa, will have support over an enor- 
mous number of vacua. Then it is impractical to solve 
Eq. ([T]) in detail. But fortunately, the probability of a 
particular vacuum is no more interesting than that of a 
specific quantum state in a macroscopic system. A more 
natural set of questions is the following: How large is the 
subset of vacua selected by cosmological dynamics, and 
how can it be characterized? 

In fact, we will ask a more restrictive set of questions. 
Let us suppose that vacua with X Xq contain no ob- 
servers, where 



A « 7.9 x 10 



-121 



(3) 



is the observed cosmological constant. Then we may re- 
strict our attention to vacua with cosmological constant 



1 In Ref. 7], A, B, ... represented specific vacua. Here they are 
indices running over all vacua. In Appendix \C\ we use indices 
a, b, . . . to run over metastable vacua only. Indices i, j, . . . are 
reserved for fluxes. 
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A of order the observed value: 



< A < An 



(4) 



(For simplicity, we exclude negative values of A; this will 
make our formulas simpler without changing our results 
qualitatively. For the same reason, we will take the upper 
bound to be sharp.) We will consider only models which 
contain a large number of such vacua, and we will be 
interested in the restriction of the probability distribution 
Pa to this subset of vacua. Thus, we ask: 

(1) How large is the set of vacua which (a) satisfy 
Eq. (j4]) and (b) are selected by cosmological dy- 
namics? 

(2) What characteristics distinguish typical vacua in 
the above set from typical vacua which merely sat- 
isfy (a)? 

Instead of solving Eq. (pj, we will approach these ques- 
tions by a Monte Carlo simulation. We select an unbiased 
sample of N initial vacua according to a generic initial 
distribution . For each initial vacuum, we simulate a 
decay chain: At every step, a new vacuum B is selected 
with a probability given by the branching ratios, t/ab, 
where B is the old vacuum. The chain will eventually 
reach a terminal vacuum, i.e., a vacuum with negative 
cosmological constant. We are interested in the penul- 
timate vacuum, which still has positive vacuum energy. 
For large enough N, the penultimate vacua obtained from 
the Monte Carlo simulation become a representative sam- 
ple of the ensemble of vacua selected from the landscape 
by cosmological dynamics. 



B. BP model and decay rates 

We will follow Ref. [l2| in modeling the string land- 
scape as a set of flux vacua, neglecting the backreaction 
of fluxes and moduli stabilization. The geometry is of the 
form Mi x X, where X is a compact 6 or 7-dimensional 
manifold, and M4 is the macroscopic four-dimensional 
spacetime. A p-brane wrapped on a p — 2 cycle in X 
acts as a membrane in M4. The number of cycles, J, and 
thus the number of membrane species, can range into the 
hundreds. Each species sources a four-form field strength 
in M4. The four-form flux is quantized, so each vacuum 
is fully specified by the set {rii}, i — 1, . . . , J, of integer 
flux numbers. The vacua form a J dimensional grid in 
flux space (Fig. [TJ. 

The cosmological constant 2 can be written as 



Here, is the charge of the i-th membrane. The sec- 
ond term is contributed by the J species of four-form 
flux. Abare subsumes all other contributions, in particu- 
lar those from vacuum loops, and will be assumed to be 
naturally large and negative: 



0(1). 



(6) 



If each flux can take several different integer values, 
then there will be an enormous number of vacua. Some 
of these vacua will accidentally have a very small cosmo- 
logical constant, even though no small scale is introduced 
into the model. If the model contains vacua with a cos- 
mological constant on the order of the observed value in 
our universe, Ao, then the cosmological constant problem 
can be solved by anthropic selection effects @, El, [H| ■ 

A more detailed model of the string landscape should 
include the backreaction of fluxes on the geometry (see 
Ref. [IH for a review). These effects are still hard to 
control, especially if we are interested in the cosmology 
of a large, representative part of the landscape. In par- 
ticular, both the charges qi and the bare cosmological 
constant would depend on the set {rii}. However, this 
is unlikely to affect a key feature of the BP model: that 
the vacua form a discretuum with tiny effective spacing 
of the cosmological constant. 

In general, a metastable vacuum can decay to any other 
vacuum. However, a typical decay will change one flux 
by one unit. Tunneling events that change more than 
one flux, or change a flux by more than one unit, can be 
regarded as a composite event. They will be much more 
suppressed than one of their constituent events and can 
thus be neglected. This simplifies our task, reducing the 
number of possible decays at each step to 2 J. The only 
question is which flux will change, and in which direction. 

The relative decay rates are given by Eq. ([2]), with 



= Ai 



-,-Bi 



(7) 



The constant Aj depends only on the host vacuum and 
drops out in the branching ratios. The instanton action 
for decreasing the magnitude of the i-th flux, m, by one 
unit is given by [l6j 



B(\m 

-Bflat(K 



M-i) 
kl - 1) 

r(x,y) 



£flat(K| -> \rii\ - l)r(x,y) 
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(5) 
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2A 
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ll; 



TY- 1 ' 



(8) 



2 We use units in which 8nG = c = h = 1. 



where A is the cosmological constant in the first vacuum. 
Upward tunneling is further suppressed by the relative 
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qjnj 




FIG. 1: A two-dimensional slice through the J-dimensional flux grid. Surfaces of constant A are J — 1 dimensional spheres. 
Each dot represents a vacuum. Starting from initial vacua with A » 1 (outermost shell), we simulate decay chains (blue lines) 
through the landscape, which terminate when a vacuum with A < is reached. Among the "penultimate" vacua with relatively 
small positive A (green/shaded region), only a small fraction is actually accessed by the decay chains, leading to a host of 
predictions. Vacua with A of order the tiny observed value lie in a much thinner shell (schematically shown as a black circle). 
A model can be ruled out if the selection effects render this shell inaccessible. 



entropy of the two vacua: 

B (\ ni \ - \m\ + 1) = (9) 
%| + l^N)exp( >+ ff iK -^) • 

For fluxes with m — it is important to take into account 
that — > 1 and 0^—1 are two distinct decay channels 
that contribute equal amounts to the total decay rate 
entering Eq. @. 



C. Quantifying selection effects 

Given the N vacua selected by our Monte Carlo sim- 
ulation, we can answer the two questions posed above. 
First, let us focus on how the selected vacua can be char- 
acterized and distinguished from an unbiased ensemble 
of vacua with small cosmological constant. 

In our toy model, the only variables we have available 
are the fluxes, n^. A list of J fluxes specifies a unique 
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vacuum, but a list of J probability distributions (one for 
each flux) defines an ensemble of vacua. Larger ensem- 
bles could be defined by throwing away more information, 
e.g., by specifying only a list of J expectation values (rij). 
However, we are interested in extracting as much infor- 
mation as possible from our simulation. Hence, we will 
consider the J discrete distributions 



Pi (n) = probability that the «-th flux, rii 
has integer value n . 



(10) 



Selection effects can be demonstrated by comparing the 
J distributions obtained from the selected vacua, with 
the J distributions obtained from an unbiased sample of 
vacua with a similar range of vacuum energy. 

Next, let us turn to estimating the total number of 
selected vacua from the small sample obtained by the 
simulation, and comparing this to the total number of 
vacua in a similar range of vacuum energy. This task is 
slightly more subtle. 

Associated with each probability distribution Pi(n) is 
a Shannon entropy, 



Si 



OO 

E 



Pi (n) log pi (n) 



(11) 



We shall treat the rii as independent random variables. 
This assumption is justified, and small corrections are 
given, in Appendix [A"l Then the total entropy is 



(12) 



i=l 



and the corresponding number of "states" (i.e., vacua) in 
the ensemble is 



N = exp S . 



(13) 



We will be interested in computing this number not only 
from the selected sample of N vacua, but also from an 
unbiased sample of vacua with a similar range of vacuum 
energy. Their fraction, 



elected 



(14) 



quantifies the degree to which cosmological selection 
thins out the landscape. 

We are especially interested in how many distinct val- 
ues of the cosmological constant, A, are contained in the 
spectrum of vacua, before and after cosmological selec- 
tion. By Eq. ([5]), changing the sign of a flux, rij — * — rij, 
leaves A invariant. To account for this degeneracy, let 
us consider a new set of distributions in which vacua are 
treated as identical if they differ only by the signs of 
fluxes: 



Pi(n) 



Pi(0) , 
Pi(-n) 



Pi{n) 



n = 0; 
n > 0. 



(15) 



From this set of pi, we can compute the number M of 
different values of the cosmological constant, by the ana- 
logues of Eqs. (P2J and (fl3|l. 

As we remarked above, we assume an initial proba- 
bility distribution that depends only on A. This ensures 
that pi(—n) = pi(n) initially (up to finite-sample effects), 
and by since the decay rates depend only on \rii\, the con- 
dition is preserved by dynamics. Hence, 



Pi(n) = Pi(-n) = 



hPi(\n\). 



n = 0: 



(16) 



For this reason, we only display Pi(n) in all figures below. 

The N vacua in our sample will be distributed over 
a much wider interval of A than that of Eq. (|H); very 
likely, none of them will have A ~ Ao. However, the 
distribution of A is smooth on scales much larger than Ao , 
so that the total number of dynamically selected values 
of A satisfying Eq. A/w, can be easily estimated: 



J\f dN 



An — 



A=0 



N d\ 



(17) 



A=0 



where the distribution dN/dX can be obtained from our 
Monte Carlo simulation after suitable binning, and, from 
Eqs. (HU and (USD, 



M = exp[-^p»(n) log£i(n)] 



(18) 



In AppendixEJwe note that the fluxes are only approx- 
imately, but not completely independent random vari- 
ables. We argue that an appropriate correction can be 
implemented by generating a new sample of N' vacua 
from the probability distributions pi(n) obtained from 
the penultimate vacua on the Monte Carlo chain, and 
replacing N by N' in Eq. (JT7J). 

In Appendix [B] we explain how to estimate the num- 
ber A7^ o nsolcctod of distinct values A - A in the model 
(before cosmological selection) . An approximate formula 
was given in Ref. [Hj but we propose a more precise 
method that better accounts for degeneracies of vacua. 



D. Initial conditions 

There is no universally accepted theory of initial condi- 
tions, and it is not our purpose to propose such a theory. 
Given initial conditions, we are interested in the effects 
of the cosmological dynamics on the selection of vacua 
with small cosmological constant. For definiteness, we 
will assume an initial probability distribution that de- 
pends only on A and favors large values of the vacuum 
energy. 

This type of initial distribution would follow, for ex- 
ample, by assuming that all vacua are equally likely. The 
number of vacua with cosmological constant less than A 
grows very rapidly, like (A — Abaro)' 7 ^ 2 - Hence, the over- 
whelming majority of vacua have large positive values of 
A. 
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FIG. 2: The number of selected vacua with A ~ Ao, log7vA Q , 
obtained from different initial ensembles starting at different 
values of (A), in Model 1. For initial values (A) > 0.8, the 
number of selected vacua stops shrinking. This indicates that 
at such high values, selection effects are negligible. Therefore, 
the details of the initial conditions are irrelevant as long as 
large values of A are preferred, as explained in the text. 
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An even stronger preference for large initial A follows 
if we adopt the tunneling proposal [IT], [18| : 



Pi 0) oc 



e- 24 - /\ A>0: 
0, A < . 



(19) 



This assigns negligible probability to vacua with small 
cosmological constant. Vacua with a large positive vac- 
uum energy are then favored not only because of their 
greater number, but because of their greater individual 
probability. 

It makes little difference which of these specific propos- 
als we adopt. The above distributions diverge at large 
values of A and must be cut off, so most initial vacua will 
be at the cutoff. We will choose a cutoff A < 1. It seems 
reasonable to ignore vacua with A > 1 since semiclassi- 
cal gravity cannot be valid in this regime. One signature 
of this breakdown is that action of instantons mediating 
their decay becomes of order unity or less. This means 
that the vacua cannot be treated as distinct, metastable 
states. 

Fortunately, our results show little sensitivity to the 
precise choice of the cutoff. We have run sets of de- 
cay chains starting from different initial values of A. As 
the initial value is increased, the selection effects at first 
become more important. This is not surprising, since 
a larger initial value implies a longer decay chain, and 
hence, more opportunities for preferred decay channels 
to imprint directional effects on the probability distribu- 
tion. Eventually, however, the distribution shows asymp- 
totic behavior, changing little as the initial A is further 
increased (see Fig. O in the model here, this occurs near 
A = 1). This is also reasonable: at large cosmological 
constant, all decay channels become unsuppressed. Di- 
rectional selection effects appear to set in only when A 
becomes sufficiently small. Therefore, the precise cutoff, 
and thus the choice of initial cosmological constant, do 
not affect our results. 



FIG. 3: The red, wiggly curves show the probability that the 
i-th flux is n, Pi(n), in the penultimate vacua of N = 7000 
decay chains starting from A > 0.8 in Model 1, as a function 
of i. For comparison, the blue, smooth curves show the prob- 
abilities that would have been obtained without cosmological 
selection, just from restricting to vacua with small cosmolog- 
ical constant. The differences are immediately apparent and 
quite drastic. For example, no fluxes greater than 4 survive 
the selection process; without selection, there would be many 
such fluxes. The fluxes associated with small charges (small 
i) are anomalously low after selection. 



This clearly alleviates the dependence of the local 
probability measure on initial conditions. However, it 
does not eliminate this dependence entirely. For example, 
if the initial probability distribution was concentrated on 
a particular vacuum (rather than dependent only on A), 
then the selected vacua would generically show a signifi- 
cant imprint of this origin. This is as it should be. If we 
reject the global viewpoint, there is no reason to expect 
that all initial conditions must be washed out. 



III. RESULTS 



We will first present results for a model with J = 250, 



A barc = -0.5, and q t = 0.01494 
The model contains A/\ unsolcctod » 



terval < A < A , with A^ nsclccted 



0.03^ ("Model l"). 3 
10 121 vacua in the in- 
10 60 distinct val- 



s Strictly, the charges should be incommensurate in order to avoid 
additional degeneracies. Since this could be implemented by ex- 
tremely small modifications of the , we ignore such degeneracies 
here. 
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ues of A. We considered initial distributions centered on 
(A) = 0.1, 0.2, . . . , 1.4. For each initial distribution (for 
example, the red/dark shell in Fig.[T]), we ran 1000 Monte 
Carlo decay chains selecting 1000 vacua with small pos- 
itive cosmological constant (in the green/light shell in 
Fig. EE]). 

The probability distribution of fluxes characterizing 
the selected vacua becomes independent of the initial dis- 
tribution for initial values of A > 0.8 (see Fig. [2]). This 
asymptotic distribution is what we are after: it is the one 
picked out by the initial conditions chosen in Sec. Ill Dl 
To improve the statistical quality of our final data, we 
therefore merged the distributions obtained from initial 
(A) = 0.8, 0.9,..., 1.4, effectively obtaining N = 7000 
runs. 

The resulting probability distribution of fluxes charac- 
terizing the selected vacua, p~i(n), is shown as a function 
of i in Fig. [21 and as a function of n in Fig.0J For compar- 
ison, we show the probabilities that would have been ob- 
tained from a random sample of vacua with comparable 
cosmological constant, without cosmological selection. 

The overall strength of the selection effects (our first 
question in Sec. Ill A[) can be quantified by the ratio of 
selected vacua to the total number of vacua with distinct 
cosmological constant in the interval < A < Ao (the 
thin black shell in Fig. [TJ. We find from Eq. (JT7J) that 



■A/x, 



10 



SO 



A/a 



yW*unselected 
Ao 



10" 



(20) 



In other words, only one out of every 10 vacua with 
small cosmological constant is actually accessed by cos- 
mological selection. 

After taking into account degeneracies, there remain 
effectively Af\ ~ 10 22 distinct values of A that are ac- 
cessed by cosmological selection. Hence, the cosmological 
constant problem can still be solved in this model, albeit 
with a much reduced discretuum density. 

In order to answer the second question posed in 
Sec. Ill Al let us inspect the probability distributions pi(n) 
in more detail, and look for characteristics that distin- 
guish the selected vacua from the full set of vacua with 
small cosmological constant. In principle, every one of 
the J fluxes will have a different probability distribution, 
leading to O(J) predictions. For example, the first flux 
n\ vanishes with 95% probability. Before selection, by 
constrast, it would have been unlikely to do so (10%; see 
the first panel in Fig. 2]). 

Instead of going through each flux, it is more illuminat- 
ing to identify general trends, paying special attention to 
features that would have been extremely unlikely in the 
unsclcctcd sample. Quite clearly, small charges tend to 
have smaller fluxes. Interestingly, this is precisely the 
opposite of what would be naively expected without cos- 
mological dynamics. 

Let us define imax(n) to be the largest i such that there 
is less than 1% chance to find one or more rij > n in 
the range 1 < j < i. For the selected vacua, we have 



W(l) = 14, w(2) = 31 and i roax (3) = 77. With- 
out cosmological selection, by contrast, it is practically 
certain that one or more fluxes in the above ranges will 
exceed n. 

Strictly speaking, we cannot resolve probabilities less 
than Af _1 , the limit from the number of Monte Carlo 
runs. However, one would expect to obtain much 
stronger predictions by involving analytic arguments. 
For example, the ^(4) graph shows that for each i < 34, 
the probability of finding rii = 4 is "zero", i.e., less than 
A~ _1 . Dynamically, this must have resulted from a strong 
preference for such fluxes to decay — a preference that 
must have been present over more than 34 steps in the 
decay chain. For example, if we assume that n\ = 4 at 
some point in the decay chain, it must be true that long 
before the end of the decay chain, the probability for it 
to remain untouched must become very small. But then 
it will be small for each of the remaining steps. Hence, 
we expect that pi(4) is not just less than A 7 ' -1 , but in 
fact exponentially small. 

We have studied a range of other model parameters 
to ensure that our results are not peculiar features of 
this model. For example, in Model 2 with J = 200, 
Abaro = -0.6, and q t = 0.01494 



0.033-, WC find 

. iTTTTry — 10~ 23 ; the selection effect is smaller but still 

considerable. The number of different values of A ~ Ao 
is _/\/ r A J o nsolcctod pa 10 39 before cosmological selection and 

A/"a rs 10 16 after selection. Thus, as in the previous 
model, cosmological selection effects yield predictions but 
keep the model viable. 

Of course, there will also be models which are ruled 
out by cosmological selection effects. The simplest way 
to find such models is to reduce the discretuum density 
before selection, for example by reducing the number J 
of fluxes or decreasing |Abare|- If the model contains only 
a few vacua with A ~ Ao to start with, chances are that 
there will be none left after cosmological selection effects 
are taken into account. For example, consider Model 3 
with qi = 0.01494 + 0.03^, J = 200 and A ba rc = -0.4. 
The discretuum density per A is A?.» nselected w 10 24 be- 
fore selection but A/a ~ 10~ 12 after selection. In other 
words, not even a single vacuum with A ~ Ao can be ac- 
cessed. This model is ruled out by cosmological selection. 

But the strength of selection effects can also vary. For 
example, Model 4 with q t = 0.01494 + 0.033- , J = 150 ' 



-1.0 has the same A/\ unsolected 

Ao 



10 24 as Model 3. 



But after cosmological selection, the discretuum density 
per Ao is Af\ — 10 13 , still larger than 1, so the model 
remains viable. 

Finally, we should mention that there can also be mod- 
els in which cosmological selection effects are unimpor- 
tant. For example, for very small charges qi <C 1, the 
typical fluxes rij will be large, and their standard devia- 
tions drii in the interval < A < Ao can be much larger 
than the typical enhancement or suppression of fluxes by 
cosmological selection effects. However, such models ap- 
pear unnatural in the context of string theory. Hence, we 
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i=l i=30 i = 50 




FIG. 4: Specific examples of probability distributions Pi(n), plotted against n for 6 of the 250 fluxes in Model 1. The 
cosmological dynamics drives most fluxes into a narrower range (red/light shaded bars) than the distribution obtained from 
landscape statistics alone (blue/dark bars). 



expect that the true string theory landscape will exhibit 
nontrivial cosmological selection effects. 



IV. DISCUSSION 

We can understand the selection behavior qualitatively 
by studying the tunneling rates, e~ B , where B(qi,m) is 
the instanton action given in Eq. ([5]). Let us neglect 
tunneling events that increase the magnitude of a flux 
rii. Then DB/dn < 0. This means that among two un- 
equal fluxes with similar charge « qf), the larger flux, 
max{ | rii |, |^j|}, is more likely to decay. The dependence 
on the size of the charge is more complicated. For most of 
the parameter space explored here, one finds dB/dq > 0, 
which implies that among two equal fluxes \m\ = \fij\, 
the flux associated with the smaller charge, mm{qi,qj}, 
is more likely to decay. (Exceptions can occur in the final 
stages of the decay process.) 

In order to understand how these effects combine and 
compete, we have plotted the lines of constant action 
in the flux-charge plane, for a number of values of the 
cosmological constant A, in Fig. [5] At the beginning of 
the decay chain, for large A, the action ranges only over 
an interval of a few, so the most unlikely decay and the 
most likely decay do not differ enormously in probabil- 
ity. In this regime, our simulation shows that even up- 
ward jumps (which are neglected in Fig. [5]) are not very 
suppressed. 

As A decreases below 0.8, a preference for the decay 
of fluxes with small qi and/or large \rii\ becomes appar- 
ent. The slope of the constant action lines indicates how 
these tendencies are balanced. Starting around A = 0.2, 
the elimination of large fluxes is strongly preferred. This 
explains the absence of fluxes | rij | > 5 in our results (see 
Fig. [31 last panel). In their absence, the last panel of 
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FIG. 5: Constant action surfaces for A = 0.8 (top- left), 
A = 0.4 (top-right), A = 0.2 (bottom-left), and A = 0.1 
(bottom-right) in the i-\rn\ plane. (Recall that the charges 
qi are ordered so as to increase monotonically with i.) The 
action B grows towards the bottom right with a line spacing 
of 1 in all four plots, corresponding to greater suppression of 
the decay. 



Fig. [5] shows a strong preference for the decay of the 
fluxes associated with the smallest few charges. At least 
qualitatively, this explains the peaks in the first few pan- 
els in Fig. |1 
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FIG. 6: The red, wiggly curves are the same as in Fig. 3. The 
blue, smooth curves show the probabilities for various flux 
values near the beginning of the decay chain, at A = 0.8. In 
the first few panels, it is apparent that small fluxes associated 
with large charges do not participate in the dynamics. This 
can be understood from the functional form of the instanton 
action. 



Throughout the decay, in all four panels of Fig. [5j the 
decay of fluxes with both small \rii\ and large qi is rel- 
atively suppressed. This suggests that such fluxes will 
be left untouched by the cosmological dynamics, and in- 
stead will be set by initial conditions. To check this con- 
clusion, let us compare the probability distribution pi(n) 
obtained from the N penultimate vacua, with the dis- 
tribution at the beginning of the decay chain. This is 
shown in Fig. [S] The first three panels (n < 2) demon- 
strate that the initial probabilities are indeed preserved 
for sufficiently large i. 



V. COMPARISON WITH AN ALTERNATIVE 
MEASURE 

We have computed the probability distribution pi(n) 
from the measure proposed in Ref. [Tj . Other cosmologi- 
cal measures have been proposed, and it is important to 
ask which, if any, of the proposals is on the right track. 
In our view, this problem should be treated no differ- 
ently from any other challenge to come up with the right 
theory describing a set of phenomena. The measure prob- 
lem, like the question of the dynamical laws governing the 
early universe, is an aspect of the theoretical challenge 
faced in quantum cosmology. It should be approached 
by a combination of seeking out simple guiding princi- 
ples (such as the causal-diamond viewpoint of the uni- 



verse that leads to the holographic measure used here), 
and elimination. If a measure, or a theory, is poorly mo- 
tivated, then it is unlikely to get much attention; and if 
a theory with a given measure yields absurd predictions, 
then the theory or the measure or both must be wrong. 

An interesting proposal which has been quite sharply 
formulated is the measure of Garriga, Schwartz-Perlov, 
Vilenkin, and Winitzki [5] (GSVW). This measure has 
been applied to a simp lified BP model by Schwartz- 
Perlov and Vilenkin [13| (SV), so let us compare results. 



A. Initial conditions and computational complexity 

Before going into the details, we note some formal dif- 
ferences. The GSVW measure depends on initial con- 
ditions only through the sign of the initial cosmologi- 
cal constant. 4 We know of no fundamental reason why 
cosmological dynamics should be required to have this 
property. It would be nice if we were spared the task 
of understanding the initial conditions. But even in the 
GSVW measure one would still need to understand why 
the universe started in a de Sitter vacuum. 

The measure we have used depends on initial condi- 
tions in additional ways, for example through the value 
of the initial cosmological constant. As we have seen, 
however, the dependence is weak if we assume that it 
enters mainly through A. But clearly, initial conditions 
remain an free variable when applying the holographic 
measure; we have argued only that within a wide range 
that appears to be physically reasonable, the measure is 
not very sensitive to them. 

It is unclear whether in practice, the GSVW measure 
can actually be applied to a rich landscape such as that 
of string theory, as we have done here for the holographic 
measure. As shown in Ref. jl3| . the GSVW probability 
distribution is dominated by the offspring of the longest- 
lived metastable vacuum. This requires finding the eigen- 
vector s with smallest-magnitude negative eigenvalue of 
a matrix with rank ten to the hundreds. Perhaps statisti- 
cal methods can be developed for this purpose, but so far 
the problem is unsolved. It is harder than the problem we 
addressed here of statistically solving Eq. ([TJ) , and it may 
be computationally intractiblc in a realistic landscape, in 
the sense of Ref. [19| ■ 

In any case, SV considered only a toy-BP-model with 
J = 7 charges and a few million vacua. Then the eigen- 
vector s can be found perturbatively in the small upward 
tunneling rates. Because of the much smaller number of 
vacua (in particular, their model does not contain any 
vacua with A ~ Ao), we cannot directly compare the re- 



4 More precisely, the measure is not defined for initial states com- 
posed exclusively of terminal vacua, and depends sensitively on 
initial conditions for initial states involving metastable A < 
vacua but no A > vacua. 
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suits of SV to our results for a landscape with ~ 10 250 
vacua. 

Even at a qualitative level, Ref. [l3| contains no predic- 
tions of likely flux distributions that could be compared 
to ours. However, there is one robust prediction of the 
GSVW measure, which is not shared by the holographic 
measure, and which appears to conflict with observation. 
We turn to this problem next. 



B. The staggering problem 

The most important feature in the results of SV is the 
"staggered" nature of the probability distribution. The 
GSVW probabilities are widely spread. At first order 
in the perturbative approximation pursued by SV, the 
logarithmic range of probabilities, logp^, appears to be 
larger than the number of vacua. (A similar conclusion 
was reached in Ref. [20] for a different landscape.) 

This behavior is in sharp contrast to our distribu- 
tion, and it is extremely problematic. As SV point out, 
if the staggered behavior persisted at all orders in the 
full string landscape, the cosmological constant problem 
could no longer be solved. Then either the landscape or 
the GSVW measure would be ruled out. 5 At first order 
in the SV approximation, staggering is generic, but one 
might hope that at higher order, a more uniform proba- 
bility distribution will emerge [201 ]. 

However, we will argue that staggering is a very gen- 
eral feature of the GSVW measure. We will use a result 
proven in Appendix [C) The GSVW measure can be per- 
fectly mimicked by the holographic measure, if the initial 
state is chosen to be the GSVW eigenvector s, which 
is dominated by the longest-lived metastable vacuum in 
the landscape, (*). This result allows us to estimate the 
Shannon entropy of the GSVW probability distribution 
by following the branching trees of Ref. Q • 

Unlike in the earlier sections, we will be estimating the 
probability of vacua directly, not as a product of proba- 
bilities of individual fluxes. The Shannon entropy is then 
given by 

Sgsvw = -^2pa log pa , (21) 
A 

and the effective number of selected vacua is A/"gsvw = 
exp(S'Gsvw)- The cosmological constant problem can 
only be solved if A/"gsvw ^ 10 120 . In a staggered prob- 
ability distribution, the most likely few vacua dominate 



5 Until wc discover alternatives to the string landscape in funda- 
mental theory, or at least until wc have exhaustively explored 
alternative measures, the more conservative conclusion would be 
to regard this problem as evidence against the GSVW measure, 
not against the string landscape. This is all the more true since 
in the holographic measure we have a well-motivated alternative 
that does not lead to this problem. 



the probability and the entropy, so Sgsvw ~ 0(1), which 
is far too small. We will now argue that this type of distri- 
bution results from the GSVW measure quite generally. 

As SV point out, the (*) vacuum, by virtue of be- 
ing long-lived, will have a relatively small cosmological 
constant, and upward tunneling will be enormously sup- 
pressed. Let us denote by e the total branching ra- 
tio for upward tunneling. By Eq. (flu]) , loge ~ — 1/A*, 
where A* is the cosmological constant of the (*) vacuum, 
so that — loge » 1. After one upward tunneling, the 
cosmological constant is larger, but still only of order 
(rii)qf ~ 0(qi/y/J), which is less than 10~ 2 in natu- 
ral models, so upward tunneling remains extremely sup- 
pressed. Following SV, let us thus use upward tunneling 
as a small expansion parameter. At zeroth order, we 
consider only downward tunneling. 

SV showed that the eigenvector s is dominated by one 
component, the (*) vacuum, and that the (*) vacuum 
will have a small enough cosmological constant to access 
only terminal vacua by a single downward decay. The 
arguments leading to this conclusion generalize to realis- 
tic models. However, to show that the GSVW entropy is 
small, we will only rely on the weaker statement that the 
GSVW-equivalent initial conditions correspond to start- 
ing from a small number of initial vacua whose downward 
decay leads to terminal vacua after very few steps. 

Let us imagine that the downward decay of any 
metastable vacuum affects each flux with equal relative 
probability, l/J. Moreover, let us pretend that decay 
chains never merge; all vacua reached by downward de- 
cays are different. These assumptions vastly overesti- 
mate the spread of probabilities and thus overestimate 
the GSVW entropy, for two reasons. First, we are at 
small cosmological constant and some decay channels will 
be much more suppressed than others; vanishing fluxes 
cannot decay downward at all. Second, the effective num- 
ber of selected vacua will also be smaller since some vacua 
will be reached by more than one decay chain. 

For simplicity we assume that the tree contains only 
metastable vacua until the D-th step, at which point all 
vacua will be terminal. Starting from the the (*) vacuum, 
D will most likely be 1, but when starting a downward 
progression after the first upward jump, D can be a few, 
so let us keep it general. 

At zeroth order, we thus have D generations of vacua. 
The d-th generation contains J d vacua with normalized 
probability pd,o ~ l/DJ d (the second index refers to the 
zeroth order). A bound on the zeroth order entropy is 
therefore 

Sgsvw.o < -JjJ2 log JTjd = 1 °g £) + — ^— lo S J > 

d=l 

(22) 

This is of order one for realistic values of J < 0(1000) 
and D ~ 0(1). 

The number of states, A/qsvw o ~ DJ~i~ , is very 
small compared to 10 120 — far too small to include vacua 
with A ~ Aq. Typically, anthropic considerations cannot 
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help in such a case, since most observers will be rare fluc- 
tuations in a vacuum with A » A . Hence, the GSVW 
measure is not capable of populating the landscape at 
leading order. 

Before turning to higher order corrections, let us pause 
to ask whether the holographic measure would have suc- 
ceeded at this stage, i.e., "at leading order", when up- 
ward tunneling is neglected. Indeed, our simulation ob- 
tains large entropy even if upward tunneling is turned 
off. However, one might be concerned that this is mainly 
the result of starting from a large ensemble of vacua with 
similar, large A, rather than just a single vacuum (natural 
though such a starting distribution may be (2ljV 

But in fact the holographic measure would produce 
large S even if we restricted P' -* to a single member 
of the initial ensemble. The only thing that matters is 
that we start from a generic initial vacuum, for which 
A is of order unity. Then there will be a large number 
D ~ 0(\/~J /q) of downward tunnclings (D ~ 0(100) in 
our model) before the cascade terminates, and Eq. (|22|) 
yields Af ~ J°(yJ/i) ; which is much larger than 10 120 in 
generic models. 

Returning to the GSVW measure, the entropy in 
Eq. (f!?2"l) receives corrections at higher order. We will 
now argue that they are suppressed by powers of e <C 1 
and so they will hardly change the zeroth order entropy. 
In particular, they will not change the conclusion that 

Agsvw « !0 120 - 

At each order, we make one upward jump and gain 
one "supergeneration" , or downward cascade of vacua, to 
which the above analysis can again be applied. Because 
upward tunneling is so highly suppressed, it will almost 
certainly originate from the largest-A vacuum in the u-th 
supergeneration and lead to a strongly preferred vacuum 
with larger A. Corrections to this assumption will be 
considered below. 

Then the GSVW entropy may again be (over- 
estimated by the normalized probability distribution 

p( u ,d) = izf^ 1 (23) 

for a vacuum in the d-th generation of the u-th super- 
generation. Thus the GSVW entropy is bounded by 6 

oo D 

Sgsvw < j^5^i>(«,d) (24) 

u=0 d=0 

D e , 1 D + l, 

~ log i + i log - + —z— log J . 

1 — e 1 — e e 2 

This differs from the zeroth order result mainly through 



Here we have included the (*) vacuum (u = d = 0) to keep the 
expression simple. Strictly, it should be excluded, but this leaves 
the result essentially unchanged quantitatively; an overall factor 
(D + 1)/D would appear on the left hand side of Eq. J25t . 



the negligible term 

elogi«l, (25) 

so the GSVW entropy remains too small even at higher 
order. 

We have neglected the w-dependence of both e <C 1 
and D. Choosing an upper bound K 1 for all upward 
tunnclings and D < O(l) for the number of downward 
generations will only overestimate the entropy but leave 
the conclusion unchanged. But at very high order, u ~ 
0(y/~J/q), e might be come O(l), D will become larger, 
and the expansion might break down. Can this change 
our conclusion? 

To see that it will not, let us stop at u of order a few, 
so that our approximations hold up to this point. Let 
us estimate the remaining contributions to the entropy 
by imagining that at this point the remaining probability 
flows equally into all vacua in the landscape, resulting in 
Af vacua with probability e u /Af. (This is almost certainly 
a large overestimate unless the landscape has no sinks or 
is extended far beyond the semiclassical regime.) But 
this contributes only of order 

e"log^ (26) 

to the GSVW entropy. In a realistic landscape, we may 
safely assume that log A/" <C e~" even for u ~ 0(1). 
Hence this contribution is also negligible. 

We conclude that the GSVW entropy is dominated by 
the zeroth order term, Eq. ()22[> . Its small size indicates 
that staggering is a general problem for the GSVW mea- 
sure. It arises because the GSVW measure, in effect, 
attempts to populate the landscape starting from the 
longest-lived metastable vacuum. This fails since it relies 
on decay channels of enormously small relative probabil- 
ity. With reasonable initial conditions (such as starting 
with a randomly chosen vacuum or an ensemble at large 
A), this problem is absent in the holographic measure. 
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APPENDIX A: STATISTICAL INDEPENDENCE 
OF FLUXES 

Let us discuss the assumption that the J fluxes are sta- 
tistically independent. We will argue that there is some 
minor interdependence because J is finite, and we will 
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explain how we correct for the corresponding constraint. 
Since the correction is small, some readers may wish to 
skip this appendix. 

If there are correlations between the distributions 
Pi(n), then their entropies do not simply add. In this 
case, the total number of states will be lower than the 
"maximum randomness" formula, Eq. (|12p . would make 
us believe. This will not be a concern for our sample of N 
initial vacua, nor for our unbiased samples of vacua with 
small cosmological constant (without including selection 
effects from cosmological dynamics). As we discuss in 
Appendix IBI below . one can generate suitable samples by 
treating the fluxes as statistically independent. However, 
it is a concern for the N vacua we arrive at by selection 
though a Monte Carlo chain. 

In fact, an obvious correlation is introduced by our 
choice to consider the penultimate vacuum in each chain. 
By construction, these vacua must have a positive cos- 
mological constant small enough that the next flux decay 
can make A negative. Thus, they satisfy the constraint 

W < 2|A baxe | + (2 nj - l)g? (Al) 



flux depends on the other fluxes only through A (and not, 
for example, though terms like (nj — rii+i) 2 , which could 
introduce clustering). The vast majority of possible flux 
combinations at each step in decay chains starting from 
a fixed value of A are mapped to a small range of A, so 
the relative decay rates determining the next step cannot 
be sensitive to the detailed distribution of fluxes. 

A very conservative upper bound on such correlations 
can be obtained by ignoring averaging. Instead, let us 
imagine that decay chains fall into two families, one of 
which has higher values of A at each step. But they can- 
not differ by more than the amount by which the final 
decay, at the end of the chain, changes A (or else, the 
steps can just be relabeled so as to match the families 
more closely). At the beginning of the chain, the decay 
is random for either family, since A is large. So the two 
families can differ at most by the effect of one additional 
decay at the end. Thus, we can bound dynamical cor- 
relation effects by comparing our sample of penultimate 
vacua with the ultimate vacua they decay to. But in the 
models we study, the corresponding probability distribu- 
tions, Pi(ri), barely differ. 



for at least one flux j. This constraint would be unim- 
portant if J was very large, but in realistic models, J is 
only of order a few hundred. 

The number of vacua with distinct vacuum energies of 
order Ao was estimated in Eq. (fT7|) . The correction due to 
the constraint (jAip can be approximately implemented 
as follows. The entropy of the unconstrained ensemble is 
larger than the true entropy, with the extra states com- 
ing from vacua that fail to satisfy Eq. (|A1|) . But consider 
the distribution of vacua over the cosmological constant, 
dM /dX, in both the constrained and unconstrained en- 
sembles. In the range of A where Eq. (|Aip is guaranteed 
to hold (which ranges from A = to an upper bound 
much larger than Ao), the two distributions dM /d\ will 
agree well. 

From this viewpoint, the problem with Eq. (|17[) is not 
that M is computed from the unconstrained ensemble; 
the problem is that dN/dX is not. But this is easy to 
correct. Let us keep only the probabilities Pi(n) obtained 
from the Monte Carlo chain, but replace the actual sam- 
ple of N vacua satisfying the constraint with a broader 
sample of N' vacua randomly generated from the p, (n) , 
with the pi(n) treated as truly independent. This will 
spread the vacua over a larger range of A and thus re- 
duce the number of vacua with A ~ Aq : 



^ = A ° -dx 



A=0 



M dN' 
Xo M~dX 



(A2) 



A=0 



In the models studied in this paper, this correction re- 
duces the number of vacua with A ~ Ao by about one 
order of magnitude. 

In principle, additional correlations could arise through 
the dynamics. By Eq. {5}, however, the decay rate of each 



APPENDIX B: GENERATING AN UNBIASED 
SAMPLE 

Consider all the vacua whose cosmological constant lies 
in a small interval SX around some value A. On a number 
of occasions in this paper, we require an unbiased sample 
of such vacua. This task is equivalent to throwing all 
the vacua in the interval SX into a bag, pulling one out, 
returning it 7 , and repeating this process N times to get 
a sample of size N. 

But the vacua are not in a bag; they are arranged 
in a grid on which the above interval corresponds to a 
spherical shell (see Fig. [I}. Because the grid has a high 
dimension J, we must be careful to avoid directional ef- 
fects that tend to pick vacua from preferred patches of 
the shell. 

To generate a sample, we first need a set of probability 
distributions for the fluxes, Pi(ri). Vacua can be gener- 
ated by picking the fluxes rij randomly according to the 
distributions pi (n) . It is not necessary that all vacua thus 
generated lie in the interval 5X. We can reject those that 
do not, and generate vacua until we have N that do. 

However, this procedure will generate directional bias 
if the ensemble of vacua defined by the pi (n) is not dis- 
tributed spherically symmetrically in flux space. For ex- 
ample, if we specified that Piin) is identical for all i, then 
the distribution would not be spherically symmetric since 
fluxes with large charge qi will tend to contribute more 



7 In the cases of interest here, the number of vacua in the interval 
<5A will be much larger than the required sample size TV, so this 
step is not important. 
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to the cosmological constant. 

A spherically symmetric ensemble can be generated as 
a kind of canonical ensemble. Let us define a "temper- 
ature", T — 1//3, dual to the cosmological constant A. 
Each flux contributes an "energy" nfqf/2, so its parti- 
tion function is 

oo , 2 2 \ 

Z 0)= J2 ^pf-^J • ( B1 ) 

iii — — oo 



Treating the fluxes as independent, we obtain a total par- 
tition function 



— /3AA 



where 



AA(ni, . . . ,nj) 



J 9 9 



(B2) 



(B3) 



Recall that A = -|Abare| + A A by Eq. (5]). Since the 
probability for each vacuum, Z~ 1 e~ l3AX , depends on the 
fluxes only through the cosmological constant, the en- 
semble is manifestly direction-independent. 

The individual partition functions, z^, can be evaluated 
analytically for 8 /3qf /2 < 1: 



(B4) 



and the expectation value of the "energy" AA is given by 

J 



(AAH/3) 



2(3 ■ 



(B5) 



In order to generate a sample containing vacua near A, 
we choose 9 



J 



2(A+|Abaro|) ' 



(B6) 



and pick the fluxes rii randomly according to the proba- 
bility distributions 



p i (n) = z j - 1 exp (-/?^) 



(B7) 



This condition is not overwhelmingly satisfied in the models stud- 
ied here, so we evaluated the partition function numerically. This 
turns out to make only a small difference, but all results in the 
main part paper were obtained in this more precise manner. We 
use the analytic approximation in the formulas for ease of pre- 
sentation. 

ft is not essential that j3 is chosen just right, since in the cases 
of interest in this paper the ensemble JB2I ) has a larger spread 
than the interval SX that we want our sample to lie in. As long 
as the distribution has significant overlap with it, a sample can 
be efficiently obtained by rejecting vacua that fail to lie in the 
desired interval. 



In this paper we have two distinct occasions to apply 
this method. The first is the task of generating initial 
vacua for our Monte Carlo chains. For a given initial A we 
generate N = 1000 initial vacua according to Eqs. (|B6[) 
and (|B7p . From them, we generate 1000 decay chains and 
inspect the final probability distribution pi(n) obtained 
from the penultimate vacua on the chains. We increase 
the initial A and repeat the process until the final distri- 
bution shows asymptotic behavior, i.e., until an increase 
in the initial cosmological constant no longer modifies our 
result. 

The second application arises in estimating how many 
vacua with < A < Ao are present in a given BP 
model, before cosmological selection is taken into ac- 
count. We will begin by discussing why the estimate 
given in Ref. [12j | is too crude in the present context. 

By Eq. ([5]), vacua with small cosmological constant 
< A < Ao lie in a thin shell of radius 



1 1/2 



and width 



Ar = A /|2A 



bare 



1/2 



(B8) 



(B9) 



in flux space. Its volume is uij^ir J ~ 1 Ar, where u)j-i = 
27r J / 2 /r( J/2) is the area of a unit J — 1 sphere. 

In the limit where r 2 3> ^2qf, directional effects in the 
lattice can be neglected, and the number of vacua in the 
shell can be estimated by comparing the volume of the 
shell to the volume of a single cell in the grid, n/=i 1i • 
In the same limit, the expectation values of the fluxes, 
(rii), will be much larger than unity. Under the stronger 
condition that r 2 ^> {^2,qi) 2 J : the probability that at 
least one flux vanishes is very small. Then most values 
of A in the interval will be 2 J -fold degenerate, and the 
number of different values of A in the interval < A < Ao 
can be estimated as 



3 /2| J / 2 



1 



Ao 



bare 



(B10) 



However, here we will study models in which the shell 
radii of interest are of the same order of magnitude as the 
diagonal of a grid cell: r 2 ~ <?? ■ F° r the same reason, 
typical vacua on the shell will have many vanishing fluxes 
(rij = 0). The A-degeneracy of typical vacua will vary but 
will be less than 2 J . Hence, Eq. (|B10|) is unreliable. 

The canonical ensemble yields a better estimate. With 
(3 chosen so that (AA) « |Ab ar e|> we obtain a canonical 
ensemble characterized by the probability distributions 
(|B7|) . In order to account for degeneracies, we treat vacua 
with the same cosmological constant as identical. The 
corresponding probability distributions pi(n), with n > 
0, can be computed from Eq. (|15|) . The full number of 
vacua in this ensemble, Jv is given by Eq. (TTS|) . The 
number of vacua in the tiny interval < A < Ao is given 
by Eq. (fTT|) . The distribution dN/dX can be obtained by 
generating N vacua randomly from Pi{n) [or Pi(n)] and 
binning them suitably. 
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APPENDIX C: GSVW-EQUIVALENT INITIAL 
CONDITIONS 

In this appendix we show that our measure reproduces 
the GSVW measure [5j if the initial probability distri- 
bution in Eq. (fT]) is taken to be a particular eigenvector 
singled out by the GSVW analysis. 

We should warn that this initial condition is a very 
unnatural choice from the viewpoint of the holographic 
measure, and our result should not be interpreted as an 
equivalence of the two measures in any natural, physical 
sense. Indeed, finding the appropriate eigenvector s is the 
most challenging task in applying the GSVW measure. 
We are in effect using a nontrivial intermediate result of 
the GSVW analysis as an initial condition for the holo- 
graphic measure. But this unified viewpoint allows us 
to reduce the difference between measures to a difference 
in initial conditions. In Sec. [Vl we use this viewpoint 
to clarify both the origin of "staggering" problem in the 
GSVW measure, and its absence in our measure. 

Let us briefly summmarize the GSVW measure. It is 
derived from the evolution of a "comoving volume frac- 
tion" , 

^ = ^2(-KAB.fB + Kba/a) , (CI) 
A 

where kab is the rate per unit time for a worldline in 
vacuum B to enter vacuum A (we neglect factors of the 
Hubble parameter). Let us restrict our attention to non- 
terminal vacua (A > 0), from whose volume fraction the 
GSVW probability distribution is ultimately computed 
(see below). We may then rewrite the above equation 
using the fact that terminal vacua do not tunnel back: 

^ = K Ab f b + ^ K bafa , (C2) 

A a 

where a, b only run over non-terminal vacua. 

In the presence of terminal vacua, the solution to 
Eq. l(C2]l is of the form 

/ a (i) = S «e- fcl ' + 4 2 ) e - fc2t + ... (C3) 

where < ki < k 2 ■ ■ ■ [5]. The GSVW probability is 
defined to be proportional to the "bubble abundance" 
produced by the slowest decaying state, the vector s^: 

PA = KAbS b 1] ( C4 ) 

b 

Since the subleading terms are irrelevant, we will drop 
the sub- and superscripts and write s = and k = k\. 

For simplicity, we first demonstrate the claimed equiv- 
alence for the probability of non-terminal vacua. Let us 
write the above equations in matrix form: 

f = (K-J)f, (C5) 



where J a b — Sab J2a KAb 1S a diagonal matrix, and K 
is the matrix with components n a b- The eigenvector s 
satisfies 

(K-J)s = -fcs. (C6) 

The unnormalized GSVW probabilities for non-terminal 
vacua are 

p = Ks . (C7) 

By Eq. the absolute decay rates are related to the 
branching ratios by 

EJ = K , (C8) 

where E is the matrix with components r/ab- Combining 
the above four equations, we obtain 

(1 - E)p = (1 - E)Ks 

= (l-E)(J-fc)s 

= (J - K - k + fcE)s = kEs . (C9) 

Comparison with Eq. (TT|), with initial condition P^ -* = 
s, reveals that the GSVW probability agrees with the 
holographic probability up to an overall factor k that 
will be absorbed by normalization: 

p oc P for P (0) = s . (CIO) 

To include terminal vacua, we extend the above argu- 
ment to vectors and matrices containing terminal com- 
ponents. A prime will denote the non-square matrix of 
decay rates from terminal into terminal vacua. First, we 
note that the GSVW probabilities are independent of any 
terminal components s' ascribed to the eigenvector: 

p=(k.s)(;)=( (j - )s ) < cii > 

Similarly, by Eq. (TT]), the holographic measure does not 
depend on the terminal components of the initial prob- 
ability distribution P^ -*. For convenience, we might as 
well set s' = 0. 

The remaining question is whether the (final) proba- 
bilities for terminal vacua agree. From Eq. ^ we have 




/ (l-E)(J-*)s \ 
_ ^ [-E'(J-fc)+K']s J 

( fcEs \ _ / E 0Us\ 
~ ^ fcE's ) ) \ ) 

= kr/s , (C13) 
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so Eq. (|C10p holds for both nonterminal and terminal 
components of the probability distribution. 

Note that our result does not contradict a crucial dif- 
ference between the GSVW measure and the holographic 
measure, namely that the former depends on absolute 
lifetimes of vacua whereas the latter knows nothing about 
them. (The former is constructed from the full matrix 
k; the latter only from its normalized version, 77.) To 



achieve equivalence, we had to impart the extra informa- 
tion contained in n to the holographic measure through a 
very special choice of initial condition, s. We emphasize 
again that from the point of view of the latter measure, 
this choice would be completely unmotivated; the equiv- 
alence is useful mainly for clarifying the differences of the 
two measures (Sec. |V]) . 



[1] R. Bousso, B. Freivogel and I.-S. Yang: Eternal inflation: 
The inside story. Phys. Rev. D 74, 103516 (2006), hep- 
th/0606114. 

[2] A. Linde, D. Linde and A. Mezhlumian: From the big 
bang theory to the theory of a stationary universe. Phys. 
Rev. D 49, 1783 (1994), gr-qc/9306035. 

[3] J. Garci'a-Bellido, A. Linde and D. Linde: Fluctuations 
of the gravitational constant in the inflationary Brans- 
Dicke cosmology. Phys. Rev. D 50, 730 (1994), astro- 
ph/9312039. 

[4] J. Garriga and A. Vilenkin: A prescription for probabili- 
ties in eternal inflation. Phys. Rev. D 64, 023507 (2001), 
gr-qc/0102090. 

[5] J. Garriga, D. Schwartz-Perlov, A. Vilenkin and 
S. Winitzki: Probabilities in the inflationary multiverse. 
JCAP 0601, 017 (2006), hep-th/0509184. 

[6] R. Easther, E. A. Lim and M. R. Martin: Counting pock- 
ets with world lines in eternal inflation (2005), astro- 
ph/0511233. 

[7] R. Bousso: Holographic probabilities in eternal inflation 

(2006), hep-th/0605263. 
[8] R. Bousso: Precision cosmology and the landscape 

(2006), hep-th/0610211. 
[9] L. Susskind, L. Thorlacius and J. Uglum: The stretched 

horizon and black hole complementarity. Phys. Rev. D 

48, 3743 (1993), hep-th/9306069. 
[10] J. Preskill: Do black holes destroy information? (1992), 

hep-th/9209058. 



[11] R. Bousso, R. Harnik, G. D. Kribs and G. Perez: Pre- 
dicting the cosmological constant from the causal entropic 
principle (2007), hep-th/0702115. 

[12] R. Bousso and J. Polchinski: Quantization of four-form 
fluxes and dynamical neutralization of the cosmological 
constant. JHEP 06, 006 (2000), hep-th/0004134. 

[13] D. Schwartz-Perlov and A. Vilenkin: Probabilities in the 
Bousso-Polchinski multiverse. JCAP 0606, 010 (2006), 
hep-th/0601162. 

[14] S. Weinberg: Anthropic bound on the cosmological con- 
stant. Phys. Rev. Lett. 59, 2607 (1987). 

[15] M. R. Douglas and S. Kachru: Flux compactification 
(2006), hep-th/0610102. 

[16] J. D. Brown and C. Teitelboim: Dynamical neutraliza- 
tion of the cosmological constant. Phys. Lett. B195, 177 
(1987). 

[17] A. Vilenkin: Boundary conditions in quantum cosmology . 

Phys. Rev. D 33, 3560 (1986). 
[18] A. Linde: Quantum creation of an inflationary universe. 

Sov. Phys. JETP 60, 211 (1984). 
[19] F. Denef and M. R. Douglas: Computational complexity 

of the landscape. I (2006), hep-th/0602072. 
[20] D. Schwartz-Perlov: Probabilities in the Arkani-Hamed- 

Dimopolous- Kachru landscape (2006), hep-th/0611237. 
[21] A. Vilenkin: Quantum cosmology and eternal inflation 

(2002), gr-qc/0204061. 



